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ABSTRACT 

In this paper we report the discovery and detailed radio/X-ray analysis of a peculiar giant radio 
galaxy (GRG) detected by INTEGRAL, IGR J14488-4008. The source has been recently 
classified as a Seyfert 1.2 galaxy at redshift 0.123; the radio data denote the source to be 
a type II Fanaroff-Riley radio galaxy, with a linear projected size exceeding 1.5 Mpc, clearly 
assigning IGR J14488-4008 to the class of GRG. In the X-rays, the source shows a remarkable 
spectrum, characterised by absorption by ionised elements, a characteristic so far found in 
only other four broad line radio galaxies. 
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1 INTRODUCTION 

Giant radio galaxies (GRGs) are the largest and most ener¬ 
getic single objects in the Universe and are of particular inter¬ 
est as extreme examples of radio source development and evolu¬ 
tion. With linea r sizes exceeding 0.7 MpcQ (e.g. lLara etal.ll200l ; 
llshwara-Chandra & Saikia! Il999u . they represent the oldest and 
long-lasting part of the radio galaxies population: assuming that the 
spectral ages of radio_galaxies are representative of their dynamical 
ages (Parma et al l 19991 ). GRGs result to have on average radiative 
ages in excess of 10 8 years. Despite this, GRGs display bright ra¬ 
dio features which are indicative of current or recent activity of the 
central engine. Understanding how they can manage to maintain an 
AGN active for so long and capable of producing radio emission 
for a long time (while the nucleus in m ost galaxies becomes qui- 
escent in a much shorter period, see e.g. lSchoenmakers et alJ[200ol 
and Parma et al. 119991) requires as much data as possible on a fairly 
extended sample of giant radio galaxies. 

An interesting new development in the study of GRGs has re¬ 
cently come from the detection of a substantial fraction of these 
objects among soft gamma-ray selected radio galaxies. Since 2002, 
the soft gamma-ray sky has been surveyed by INTEGRAL/ISGRI 
and Swift/BAT at energies higher than 20 keV; up to now, various 
all sky catalogues have been released (see for example. iBird et ah! 


1 Scaled for a cosmology with Hq= 71 km s 1 Mpc -1 , £2 m =0.27, £1a=0.73; 
this value is scaled from the original l Mpc to compensate for recent 
changes in the cosmological constants. 


120101 ; iBaumgartner et al]l2013l) revealing a large population of ac¬ 
tive galaxies. Around 1 % of these soft gamma-ray selected AGN 
display a double-lobed morphology typical of a radio galaxy. What 
is intriguing and particularly interesting is that a large fraction 
(—30%) of these radio galaxies have giant radio structures (Bassani 
et al. in preparation). 

The fraction of GRGs in radio sour c e sam ples is gener¬ 
ally small (~6% in the 3CR catalogue, llshwara-Chandra & Saikial 
119991) . and it becomes relevant among radio loud soft gamma-ray 
selected objects. Therefore it suggests a link between the central 
AGN properties and the extended radio structures. To date, about 20 
objects within the INTEGRAL and Swift surveys have been identi¬ 
fied as GR Gs; a few of these are new discoveries like IGR J17488- 
2338 dMolina et alj|2014l) and IGR J14488-4008, the source dis¬ 
cussed in the present paper. These findings confirm that observa¬ 
tions at soft gamma-rays are useful tools to find radio galaxies with 
giant structures but also suggests that combining radio observations 
with measurements at high energies is a new and useful approach 
to study GRG. 

Here we report on the discovery and consequent radio and X- 
ray analysis of the INTEGRAL source IGR J14488-4008, recently 
classified as a Seyfert 1.2 galaxy at redshift z=0.123. The radio 
data indicate that the source has a Fanaroff-Riley FRII morphol¬ 
ogy dFanaro ff & Rilevi fT974l ). and linear size exceeding 1.5 Mpc, 
making it a newly discovered giant radio galaxy. In the X-ray, the 
source shows a peculiar spectrum, characterised by ionised absorp¬ 
tion, a characteristic which has so far been found only in four other 
broad line radio galaxies. 
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2 HIGH ENERGY DISCOVERY AND IDENTIFICATION 

IGR J14488-4008 was first reported as a har d X-ray emitter in 
the 4 th INTEGRAL IBIS catalogue jBird et alJl20l(j) and located, 
with a positional uncertainty of 4.5 arcmin, at RA(J2000) = 14h 
18m 50.2s and Dec(J2000) = -40d 08m 31s. The source is re¬ 
ported in the Swift-BXT 70-month catalogue and named SWIFT 
J1448.7-4009. In order to reduce the positional uncertainty associ¬ 
ated with the INTEGRAL detection and therefore to favour iden¬ 
tification, in 2011 the source was the su bject o f a fol low-up cam¬ 
paign carried on with Swift-X RT jMai izia et al.ll201 lh . Within the 
ISGRI/BAT error circles, XRT detected two X-ray sources (see 
Figure 0 }: 1) a bright source (Nl), located at RA(J2000) = 14h 
48m 50.98s and Dec.(J2000) = -40d 08m 47.12s (with a positional 
uncertainty of 3. 83 arcsec) and with an X-ray flux of 2.4x l0~ 12 
ergs -1 cm -2 (see lMalizia et alj|201 ll and lMolina et alj|2012l) : 2) a 
dimmer object (N2) located at RA(J2000) = 14h 48m 50.82s and 
Dec(J2000) = -40d 10m 56.70s, with a very soft spectrum (unde¬ 
tected above 3keV). As evident in Figure Q] which is a SUMSS 
(Sydney University Molonglo Sky Survey; fBock et al .1 19991) image 
of the sky region around IGR J14488-4008, both objects appear 
overlaid to a bright and complex radio source, showing a nuclear 
region plus a two-sided lobe structure. As reported bv IMolina et al .1 
d2012h . this radio source was observed at several radiojxequencies: 
at 1.4 GHz with a flux of 48mJy (NVSS: ICondon et al.|[l998h : at 
4.85 GHz with a flux of 48 mJy _(Combined Radio All-sky Targeted 
Eight GHz Survev-CRATES: fHealev elafll20071) ; at 8.4GHz with 
a flux of 45.9 GHz and at 8.6 GHz with a flux of 46mJy (ATPM, 
iMcConnell et alJl2012t) . IGR J14488-4008 was also observed dur- 
ing the course of an optical follow-up programme by Mas etti et al.l 
d2013h : these authors found source Nl to be a Seyfert 1.2 galaxy at 
z = 0.123, with a black hole mass of (3.8±1.9)xl0 s M o and source 
N2 to be a G-type star with no peculiarities and unlikely to produce 
high energy emission. 

As can be seen from Figure |T] the resolution of the 843 MHz 
SUMSS image is not adequate to clearly identify the various com¬ 
ponents of the radio structures that characterise IGR J14488-4008 
in this band and also is unable to associate source Nl to the core 
of the radio object. IGR J14488-4008 belongs to a sample of giant 
radio galaxies detected by INTEGRAL which is currently under 
study with the Giant Metrewave Radio Telescope (GMRT) (see[3j. 
To further study the source in the X-ray band, we also carried out 
an XMM -Newton observation in AO 12. 


3 RADIO OBSERVATIONS AND IMAGING 

The GMRT observations were performed at 610 MHz (on 9 June 
2014) and at 325 MHz (on 14 May 2014), for a total of 5 hours 
at each frequency. The observing bandwidth at each frequency is 
64 MHz, split into 256 channels with Av=250kHz. The starting 
steps of the data reduction, i.e. RFI removal and bandpa ss cali¬ 
bration, were done using flagcal iPrasad & Chengalut)l2012l) . while 
self-calibration and imaging were carried out using the NRAO As¬ 
tronomical Image Processing System (AIPS). We performed self¬ 
calibration using all the sources in the GMRT field of view at each 
frequencies. At both frequencies we produced final images at full 
and low resolution, to properly image both compact and extended 
features of the radio galaxy. The noise in the final images is in the 
range lcr ~0.12-0.20 mJy b -1 at 325 MHz and lcr ~ 40 - 80/rJy 
b -1 at 610 MHz, depending on the resolution and on the position in 
field. 



Figure 1 . 843 MHz SUMSS radio image of the region containing IGR 
.114488 4008. The larger, green circle marks the BAT position and relative 
uncertainty (5.07 arcmin); the red circle marks instead the ISGRI position 
and its uncertainty (4.5 arcmin). The two objects marked as Nl and N2 are 
the two X-ray sources detected by XRT. N1 is the true counterpart to the 
high energy source. 


Figure [2] shows the 610 MHz GMRT contours overlaid on the 
same frequency GMRT image. The high resolution and sensitivity 
of the GMRT images reveal important details on the radio emission 
of IGR J14488-4008, compared to the image provided by SUMSS. 
The central component in SUMSS is actually a blend of three com¬ 
pact components. The overlay of the radio image with the XMM- 
Newton image (see Figure [3} clearly shows that the core of the 
radio galaxy is the southernmost one, and it is coincident with the 
X-ray source labelled Nl in Figure |T| None of the other two com¬ 
pact components is optically identified on the red plate of the red 
Digitized Sky Survey DSS2, and the available data (in the optical 
and in the radio) do not suggest that these two components are the 
lobes or hot spots of a FRII background radio source. Based on our 
GMRT images their spectral index or®® 2 JJjfc is 0.5, typical of com¬ 
pact radio sources. Moreover, the eastern extension of the northern 
lobe, as visible on SUMSS, is actually a background radio galaxy 
with FR II morphology (no redshift is available in the literature for 
the optical counterpart). 

A noticeable feature of the large scale structure of IGR 
114488-^1008 is the misalignment in the morphology of the lobes. 
In particular, the high surface brightness regions of the lobes (in¬ 
cluding the hot spots) are symmetric with respect to the core, and 
are aligned in position angle of —12°, while the low surface bright¬ 
ness parts are symmetric with respect from the core but their align¬ 
ment is —5°. This is suggestive of a X-shaped radio galaxy. If con¬ 
firmed, IGR J14488-4008 would be the second giant X-shaped ra¬ 
dio galaxy iSaripalli et al.l200§) . 

Table Q] reports the flux density of each component at the 
325 MHz and 610 MHz, and the spectral index between these two 
frequencies. The total radio power is typical for FR II radio galax¬ 
ies, LogP 3 25 mhz=25.36 W Hz - 1 , and the largest linear size of 
IGR J14488-4008, estimated as the projected distance between the 
peaks of the hot spots is 1.54 Mpc. This size fully qualifies IGR 
J14488-4008 as a giant radio galaxy. The core is inverted in this 
frequency range, and the integrated spectrum of the source is flatter 
than the canonical values 0.7-0.8 (assuming S v ocy - “) usually re- 
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Table 1. The radio images 


Component 

S325 MHz 

mJy 

S610 MHz 
mJy 

>10 MHz / a y 
“325 MHz w 

Core 

20.7± 0.8 

32.3± 1.3 

-0.71 ±0.14 

Lobe N 

201.9±10.1 

148.1 ± 7.7 

0.49±0.16 

Lobe S 

354.4±17.7 

263.0±13.2 

0.47±0.16 

Total 

576.9±28.8 

443.4±22.2 

0.42±0.15 


Notes: (a) The sign convention is Soc v a . 
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Figure 2. Radio 610MHz GMRT full resolution image of IGR J14488- 
4008 with radio contours (at the same frequency) overlaid; the restoring 
beam is 8.77"x4.92" with a rms of 30mJy/beam around the source edges 
and of 50 mJy/beam close to the source. 


ported for FRII radio galaxies. This is most likely due to the fact 
that the flux density of the radio galaxy is dominated by the hot 
spots and the beginning of the backflow of the lobes. We point out 
that the frequency range of the observations presented here is be¬ 
low 1 GHz, where the spectrum of the radio galaxy may still be 
dominated by the initial population of relativisti electrons, prior to 
ageing. We estimated the intrinsic physical parameters of the ra¬ 
dio lobes under equipartition assumptions (filling factor ®=1 and 
energy ratio between protons and electrons k = 1), and the values 
are reported in Table [2] The two lobes are very similar, and the 
equipartition parameters are considerably lower (at least one order 
of magnitude) compared to radio galaxies of smaller siz^l, and are 
lower than found in the two other giant radio galaxies 3C 35 and 
3C 223 under similar assumptions dOrrii et al]20ld) . 

The lack of information on the spectral index shape and high- 
frequency cutoff of the two lobes in IGR J1448-4008 does not al¬ 
low an estimate of the radiative age of the source. 


2 Such values have been derived under the usual assumption of identical 
energy in electrons and protons, and a volume homogeneously filled with 
relativistic plasma. 


Table 2. Physical parameters 


Component 

u min 

P • 

1 min 

B e q 

Lobe N 

4.5xl0~ 14 erg cm -3 

2.4xl0' 14 dyncm~ 2 

0.7//G 

Lobe S 

5.5xl0~ 14 erg cm' 3 

2.9xl0' 14 dyncm' 2 

0.8//G 


4 X-RAY TIMING AND SPECTRAL ANALYSIS 

IGR J 14488^1008 was also recently observed by XMM-Newton on 
Feb. 16, 2014, during revolution 2599. The EPIC-pn llTurner et al.1 
l200ll) data were reprocessed using the XMM-Newton Standard 
Analysis Software (SAS) version 12.0.1 and employing the lat¬ 
est available calibration files. Only patterns corresponding to sin¬ 
gle and double events (PATTERN<4) were taken into account; the 
standard selection filter FLAG=0 was applied. The EPIC-pn nom¬ 
inal exposure of 28 ks was filtered for periods of high background, 
resulting in a cleaned exposure of —17.6 ks. Source counts were ex¬ 
tracted from a circular region of 29.5" radius centred on the source, 
while background spectra were extracted from two circular regions 
of 20" of radius close to the source. The ancillary response ma¬ 
trix (ARF) and the detector response matrix (RMF) were generated 
using the XMM-SAS tasks arfgen and rmfgen and spectral chan¬ 
nels were rebinned in order to achieve a minimum of 20 counts per 
bin. No pile-up was detected for this source. The acquisition of the 
XMM data allowed us to put firm constraints on the source iden¬ 
tification with the bight radio galaxy; it also allowed us to exclude 
X-ray emission from either the lobes or hot spots of the GRG. 



Figure 3. 0.3-10 keV XMM-Newton image of IGR 114488^1008, with the 
610 MHz radio contours overlaid as in figure[2] 
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As a preliminary step in the analysis of IGR J14488-4008. 
we focused on the issue of variability, both on short time scales, 
through the analysis of the XMM light curves, and on longer time 
scales, employing light curves from ISGRI/ISGRI and Swift/BAT. 

The XMM light curves were produced in three different en¬ 
ergy bands (soft, 0.3-2 keV; hard, 2-10 keV; total, 0.3-10 keV) and 
then analysed with the XR0N0S package. The statistical analysis of 
the light curves in the chosen energy bands, does not point to any 
significant evidence of short term variability, as shown by the rel¬ 
ative x 2 - In feet, for a 100s binned light curve, we find for the soft 
band (0.3-2 keV) light curve ay 2 /d.o.f of 112.7/165, for the hard 
band light curve (2-10keV), we find a^ 2 of 121.9/163 and finally 
for the 0.3-10 keV light curve we find x 1 - 120.9/164. For the three 
light curves we find a probability of constancy of more than 99%. 
Also the analysis of the soft vs. hard light curve hardness ratio ex¬ 
cludes any kind of spectral variability. 

We checked as well the hard X-ray light curves obtained by 
IBIS/TSGRI and Swift/BAT to examine if long-term variability is 
instead present; also in this case, we found that variability on long 
time-scales (months) is not significant and could be accounted for 
by introducing a cross-calibration constant between instruments 
when fitting the broad-band spectrum, as discussed in section l4~2l 

Spectral analysis was performed using XSPEC v. 12.8.0 
(Arnaud H~1996l) ; errors are quoted at 90% confidence level for one 
parameter of interest (Ay 2 =2.71). In the following analysis, to¬ 
gether with the proprietary XMM observation, we use data col¬ 
lected in the fourth IBIS survey feird et alj l2010h . which con¬ 
sists of all exposures from the beginning of the mission (Novem¬ 
ber 2002) up to April 2008 (revolutions 12 through 530). The 
IBIS/TSGRI data analysis and the average source spectra extrac¬ 
tion has been obtained following the procedure described by 
iMolina et alj d2013l) . To enhance the available statistics, above all 
at high energies, we alsojnake use of the accessible Swift/BAT 
spectrum dBaumeartner et al .120131) of IGR J14488-4008, obtained 
from the latest 70-month cataloguqj. XMM data have been fitted in 
the 0.3-10 keV energy range, while the ISGRI and BAT spectra 
have been fitted in the 20-110keV and 14-110keV bands respec¬ 
tively. Soft and hard X-ray data have then been combined together 
in order to obtain a broad-band coverage in the 0.3-110 keV energy 
range. 

4.1 The soft 0.3-10 keV XMM-Newton spectrum 

The EPIC-pn data were initially fitted in the restricted 2-10 keV 
energy range, in order to first characterise the primary continuum. 
As a first step, we fitted the data with a simple power law ab¬ 
sorbed only by the Galactic column column density (wa g *zpo in 
XSPEC, with N® al =0.07xl0 22 cnT 2 ); as expected, the fit is very poor 
(y 2 =457.90 for 379 d.o.f.) as some additional features (such as a 
cold absorption and an iron line component) could be present. We 
therefore added first a cold absorption component to our model 
(wa g *zwa*zpo in XSPEC); the fit improves (y 2 =403.92 for 378 
d.o.f.) and the component is highly required by the data (at more 
than 99% confidence level). We then included a Gaussian com¬ 
ponent to the fit (wa g *zwa*(zpo+zga)): the fit further improves 
(y 2 =380.20, 375 d.o.f.) and the new component is required at more 
than 99% coincidence level. The iron line is indeed well detected 
at 6.45 keV and likely to be broad (we have an upper limit on the 
line width at 139 eV and an equivalent width of 102 eV). After 

3 http://swift.gsfc.nasa.gov/results/bs70mon/ 


characterising the spectrum in the 2-10 keV energy range, we ex¬ 
tended our analysis to the softer energies. As shown in Figure [4] 
evidence of a strong soft excess is visible below 2 keV. As a first 
step towards the characterisation of this soft excess, we analysed 
the data using standard thermal models. When a blackbody com¬ 
ponent (wa g *zwa*(zbb+zpo+zga)) is employed to fit our data, 
the resulting ;^ 2 is quite poor (1101.96 for 588 d.o.f.) and the pa¬ 
rameter values are inconsistent with those expected for AGN. We 
subsequently substituted the blackbody component with the XSPEC 
model mekal, but, although the fit is better (y 2 = 699.27, 588 d.o.f.), 
the spectrum still shows residuals at soft energies, where absorption 
edges might be present. We also tested if cool Comptonisation or 
reflection models could properly fit the data, but in both cases we 
did not obtain acceptable fits. In the first case, we tried to fit the soft 
excess by employing the compTT model, while the continuum was 
modelled with a power law plus a reflection component due to the 
accretion disc (local model reflionx), intrinsic cold absorption 
and a Gaussian component to fit the Fe Kcr line. The fit was how¬ 
ever poor, yielding ay" 2 of 2468.62 (for 586 d.o.f.). We then substi¬ 
tuted the cool Comptonisation model with a reflection model, con¬ 
sisting in a relativistically smeared reflector (the final model was 
wa g *(kdblur(zpo+reflionx)+zga) in XSPEC terminology), but 
again the fit was poor, having a^ 2 of 2037.42 (588 d.o.f.). There¬ 
fore, we investigated on the possibility that the observed soft excess 
could be an artefact of smeared absorption from partially ionised 
material as described bv lGierlinski & Doneil2004 ). To test this sce¬ 
nario, we have fitted the EPIC-pn data with the smeared absorption 
from a partially ionised wind from an accretion disk (swindl in 
XSPEC) in addition to the absorbed power-law model for the con¬ 
tinuum and the Gaussian line for the Fe kcr line. The parameters 
of the smeared wind model are the absorbing column density, the 
ionisation parameter (£ = L/nr 2 J3 and the Gaussian sigma for ve¬ 
locity smearing in units of v/c. Although it has been shown that 
this model requires non physical, unusuall y high terminal veloci¬ 
ties of the outflow dSchurch & Do ne 2008), it can still be used as 
a proxy for multiple-absorber, partially-covering, ionised absorber 
models that can account for the soft excess. The smeared absorp¬ 
tion model provided a good fit to the spectrum, yielding ax 2 of 
759.48 for 587 d.o.f. and the new component is highly required by 
the data (at more than 99% confidence level). In the smeared wind 
model framework, lines from atomic recombination are expected 
to be found in the source spectrum. Indeed, although the fit is ac¬ 
ceptable, inspection of the residuals show the presence of absorbing 
features at around 0.7 and 0.9 keV. Hence, we tested our data for the 
presence of the Ovu and Fe X x absorption edges, located at around 
0.73 keV and 0.96 keV respectively. The Ovu absorption edge at 
0.73 keV has r max of 2.58 /q^, while the Fexx absorption edge at 
0.96keV has T max =0.71/g 3l j’; both components are required by the 
data at more than 99% confidence level. The 2-10 keV flux and lu¬ 
minosity are 5.26xl0~ 12 ergcnr 2 s” 1 and 1.93xl0 44 erg s _1 respec¬ 
tively. To conclude, we find that the 0.3-10 keV best-fit model to 
the XMM data is wa g *zwa*zedge*zedge*swindl* (zpo+zga). 

4.2 Broad-Band Spectral Analysis 

Extrapolation of the XMM best fit model to the 20-110keV band 
provides a slightly higher flux than those obtained by fitting the 
ISGRI and BAT spectra individually. This suggests that long term 

4 In the formula L is the source luminosity, n is the density and r the dis¬ 
tance between the ionising source and the absorbing gas 
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Figure 4. 0.3-10 keV EPIC-pn spectrum and residuals,; the model em¬ 
ployed is wa g *zwa*(zpo+zga). As it is evident, a prominent soft excess 
component is present below 1 keV. 


Table 3. 0.3-110 keV Best fit spectral parameters. The model employed is 
const*wa g *zwa*zedge*zedge*swindl*(pexrav+zga)) 


Cold Neutral Absorption 
Ovu Abs. Edge 

Fexx Abs. Edge 

SWIND1 Param. 


Continuum slope 
Cut-off Energy 
Reflection 
Iron Line 


Cross. Calib. Const. 


N H =(0.17+“$)xl0 22 cm- 2 
E=0.73 keV (fixed) 

_ _ o C9+O.41 

‘max - —0 41 

E=0.96 keV (fixed) 

r _0 71 + 0.36 
1 max - u. / 1 _oji 

N H =(7.85+!*)xl0 22 cm- 2 
Log©=2.44+“;® erg cm" 1 s _1 

rr-0 ?2 +013 
^ - 0.08 

r=i.7i +0 ^ 

E c =67jJjj^keV 
R= 1 (fixed) 

E=6.45!°;“keV 
cr=0.03 keV (fixed) 
EW=93+” eV 


Cisgri=0.7U°; 4 


flux variability could be present, although one must consider the 
possibility of a more complex spectrum when high energy data are 
considered. The long term variability can be accounted for using 
a cross-calibration constant between instruments when fitting the 
broad-band (0.3-110 keV) data. When this is done using the XMM 
best-fit model, the fit is acceptable yielding a )f of 622.68 (604 
d.o.f.) and the cross-calibration constants between XMM/BAT and 
XMM/ISGRI are 0.75 ±0.21 and 0.55±0.24 respectively. However, 
inspection of the residuals indicates a curvature in the spectrum 
above 20keV, suggesting that some additional components might 
be required at higher energies to better represent the source spec¬ 
trum. The presence of a high energy cut-off was consequently tested 
by substituting the simple power law in the baseline model with 
an exponentially cut-off power law reflected from neutral material 
(wa g *zwa*zedge*zedge*(pexrav+zga) in XSPEC terminology). 
The reflection component was added to account for the presence 
of the observed characteristics of the iron line; due to the statis¬ 
tical quality of the high energy data, we fixed the reflection frac¬ 
tion to 1 given the value of the Fe line EW (since EW/R~100- 
130eV; IPerola et aD 120021 .) The fit is acceptable, yielding a x 2 
of 772.76 for 605 d.o.f. the power law slope is slightly steeper 
than the one found for the soft X-ray fit alone, being T=1.76^'®, 
and the high energy cut-off is unfortunately not very well con¬ 
strained (E cilt =75^j 8 keV). The Ovu and Fexx absorption edges 
energies have been fixed at the values found for the soft X-ray fit, 
at 0.73 keV and 0.96 keV respectively, while their depth has been 
left as a free parameter. The Ovu absorption edge at 0.73 keV has 
r m ax of 2.58^'^, while the Fexx absorption edge at 0.96 keV has 


r =0 71 +0 - 36 
‘ max 1 L —0.31 * 


The fit we obtain is our best fit, with a^ 2 of 614.99 for 603 
d.o.f (A^ 2 = 1.02). The main spectral parameters are summarised in 
Table[3] while Figure [5] shows the fitted spectrum of IGR J14488- 
4008; note that in this case the cross calibration constants are closer 
to unity as listed in Table 3). 

The source 20-100 keV luminosity is ~3xl0 44 erg/s. G iven 
the source black hole mass (Mbh=3.8x10 s M q ; Masetti et_al j2013l) 
and employing the relation used by iMushotzkv et al. feOOa )and 
scaled to the 20-100 keV ISGRI luminosity (L Bo i=25xL 2 o-iookcv). 
we estimate an average bolometric luminosity of ~6xl0 45 erg/s 
and consequently an average Eddington ratio of 0.16, making IGR 
114488^1008 a quite efficient accretor; the uncertainty on the black 
hole mass gives a small range (0.1-0.3) for the Eddington ratio 
which therefore reinforces our conclusion that IGR 114488-4008 
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Figure 5. 0.3-110 keV XMM/BAT/ISGRI best-fit plot. Black data points 
are the XMM measurements, green data points are the BAT spectrum and 
red data points are relative to the ISGRI one. 


is an efficiently accreting source. Even using the 2-10 keV lumi¬ 
nosity and adopting the bolometric correction for radio loud AGN 
discussed bv lRunnoe et al .1 l l2012l) the Eddington ratio remains high 
( 0 . 1 ). 


5 DISCUSSION AND CONCLUSIONS 

IGR 114488-^1008 belongs to a sample of radio galaxies detected 
by INTEGRAL and currently under study with the Giant Metre- 
wave Radio Telescope (GMRT). The high resolution provided by 
the GMRT has allowed a proper determination of the nature of IGR 
114488-^1008, confirming that it is a newly discovered giant radio 
galaxy (GRG) with a linear size exceeding 1.5 Mpc and a possi¬ 
ble X-shaped morphology. Within the sample of radio galaxies de¬ 
tected by INTEGRAL, IGR J14488-4008 is the second source to 
be class ified as a new GRG (the other being IGR J17488-2338 
iMolina et alj l2014h ; this confirms that the hard X-ray selection 
is indeed a valuable tool to identify this peculiar class of radio 
galaxies, albeit this might be the result of observing in this energy 
band the most powerful objects. Both IGR 114488-4008 and IGR 
J17488-2338 appear to be characterised by high accretion rates, 
but when comparing their X-ray characteristics, these sources are 
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quite different: while IGR J17488-2338 does not have a very com¬ 
plex spectrum (a part from a cold absorber partially covering the 
nucleus), IGR J14488-4008 is characterised by the presence of a 
strong warm absorber, which is indeed a very peculiar feature. 

At soft energies, more than half of the radio quiet Seyfert 1 
galaxies show complex intrinsic absorption and/or emission lines 
associated with the presence of photoionised gas, the so-c alled 
warm absorber (e.g. lBlustin et alj2005l : lMcKeman et al.l2007i) . Up 
to very recently, their radio loud counterparts, i.e. broad line radio 
galaxies (BLRG), were not expected to show signs of ionised ab¬ 
sorption. This was ascribed to the presence of a jet seen at small 
angles with respect to the line of sight that would mask the sig¬ 
nature of ionised gas with its Doppler-booster, non-thermal radi¬ 
ation. However, with the advent of more sensitive X-ray observa¬ 
tories, this assumption has been proven to be incorrect and there 
is now a handful of sources for which the presence of a warm 
absorber has indeed been confirme d (e.g. 3C 382, 3C 390.3, 3C 
445, iTorresi et alJl2012l: 4C 74.26, lBallantvndl200l . Wi t hin this 
framework, IGR J14488-4008 can be considered as a new addition 
to this small but growing sample of BLRG with complex absorp¬ 
tion signatures in their X-ray spectra. Our analysis has also shown 
that the warm absorber observed in IGR J14488-4008 has similar 
characteristics to the warm absorbers seen in other BLR G: sam e 
column densities and ionisation parameter (see lTorresi et al ,l2012h . 
although the gas outflow velocities are slightly higher, probably due 
to the different model adopted here. 

From the available information it is possible to evaluate the 
location of the warm gas in IGR J14488-4008, warning the reader 
of the many uncertainties which effect the performed calculation 
(imperfect modelling, large errors on some parame ters, v arious 
adopte d assumptions). Following the reasoning of ITorresi et al.l 
i hoii) . we estimate the minimum and maximum distance of the 
warm absorber to be R min =(2GM)/v 2 and Rmax = Lj on /(£N H ). Since 
the ionising luminosity (L; on ) derived in the range 1-1000 Ryd¬ 
berg (13.6eV-13.6keV) from our best fit is 2.9xl0 44 erg/s, while 
the other parameters required in the above formula are those re¬ 
ported in Table 2, we estimate the location of the warm gas to 
be between 0.01 and 4.4 parsecs. In IGR J14488-4008 the loca¬ 
tion of the Broad Line Region (BLR) and the torus can be ob¬ 
tained following iGhisellini & Tavecchiol <120081) and ITorresi et al.l 
j2012l) . Assuming Li on =Ldi s k, we find RBLR=10 l7 L^ s 2 k =0.017 pc and 
Rtorus=2.5xl0 18 L* i /2 k =0.4pc, where L disk is expressed in units of 
10 45 erg/s. These values suggest that the bulk of the highly ionised 
absorbing gas detected in IGR 114488-^1008 is located internally 
to the torus. We also note that among BLRG with a warm ab¬ 
sorber component, only for 3C 445 a similar location is reported 
ITorresi et al.ll2012h . while in other objects the warm gas is found 
beyond the torus. IGR 114488^1008 is also similar to 3C 445 for 
the high outflow mass. 

In summary, IGR J 14488^-008 is a remarkable AGN; firstly 
discovered to be a high energy emitter by INTEGRAL, it has been 
identified as a Seyfert 1.2 galaxy thanks to an optical follow-up 
study and then further characterised as a Giant Radio Galaxy by 
means of a dedicated programme which is being carried out at the 
GMRT facility. IGR J 14488^1008 appears to be a highly efficient 
accretor, and is also peculiar in the X-ray domain, being one of the 
few broad line radio galaxies known to have a warm and highly 
ionising gas in its nuclear regions. Unlike other BLRG with warm 
absorbers, the warm gas in IGR J14488-4008 is located well within 
the torus region, making it similar to only another BLRG, 3C 445. 
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